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Abstract 



We study the 3-point functions generated at recombination in the squeezed triangle limit, when one 
mode has a wavelength much larger than the other two and is outside the horizon. The presence 
of the long wavelength mode cannot change the physics inside the horizon but modifies how a 



i-^ ■ late time observer sees the anisotropies. The effect of the long wavelength mode can be divided 
^p" 1 . into a redefinition of time and spatial scales, a Shapiro time delay and gravitational lensing. The 
separation is gauge dependent but helps develop intuition. We show that the resulting 3-point 
function corresponds to an /nl < 1 and that its shape is different from that created by the /nl (or 
■ local) model. 

> 

1 Introduction 

Conventional slow-roll inflation predicts that the primordial perturbations are Gaussian to an ex- 
treme level of accuracy, with deviations of order 10 -6 0121- Any departure from the simplest picture, 
like the addition of light scalars during inflation or models with modified inflaton dynamics, leads to 
an enhancement of the non-Gaussianity and many of these models predict a 3-point function which 
should be detected by the forthcoming experiments [31 QJ El EH El • 

The relation between the primordial fluctuations and what we observe today is in general non- 
linear. Although the linear approximation is quite good for studying the CMB anisotropies on 
large enough scales, small non-linearities will generate higher order moments even starting from a 
completely Gaussian primordial spectrum. The leading sources of non-Gaussianity are related to 
secondary anisotropies, anisotropies imprinted after recombination {z ~ 1100). These anisotropies 
can be divided into scattering secondaries, when the CMB photons scatter with electrons along 
the line of sight, and gravitational secondaries when effects are mediated by gravity (for a review 
of all the different secondary effects see [HI). Examples of scattering secondaries are the thermal 
Sunyaev-Zeldovich effect, where hot electrons in clusters and other collapsed objects transfer energy 
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to the CMB photons, the kinetic Sunyaev-Zeldovich effect produced by the bulk motion of the 
electrons in those clusters, the Ostriker-Vishniac effect, produced by bulk motions modulated by 
linear density perturbations, and effects due to the patchy nature of the reionization process. The 
various correlations between these different effects lead to non-zero 3-point functions and other higher 
order effects (see for example JJ}_). As the growth of structure proceeds, density inhomogeneities, 
bulk and thermal motions grow and become quite large on small length scales. As a result the 
scattering secondaries are produced mostly at late times and are most significant on small angular 
scales. 

Gravitational secondaries arise from two separate effects, the change in energy of photons when 
the gravitational potential is time-dependent and gravitational lensing. Gravitational secondaries 
are different from scattering secondaries in that even as structure formation evolves the depth of 
gravitational potential wells does not grow significantly, always being of order 10~ 5 ( 1 ). At late times, 
when the Universe becomes dominated by the cosmological constant the gravitational potential on 
linear scales starts to decay. This leads to the largest of the gravitational secondaries, the integrated 
Sachs- Wolfe effect (ISW). The ISW effect is linear in the gravitational potential and affects mainly 
large angular scales. Other secondaries that result from a time dependent potential are the Rees- 
Sciama effect, produced by the time evolution of the potential on non- linear scales, and the moving 
cluster effect, where time-dependence is produced by motions. Both these effects are subdominant. 

The fact that the potential never grows appreciably means that most second order effects created 
by gravitational secondaries are subdominant to those created by scattering ones. The exception are 
effects related to gravitational lensing. The reason is that the total defection angle of a photon on 
its way from the last scattering surface can be more than an order of magnitude larger than <j). This 
is because for a mode of wavenumber k, there are k/H independent regions along the line of sight. 
Gravitational lensing conserves surface brightness so it does not create anisotropics, it only modifies 
existing ones. Lensing is thus a second order effect but the resulting anisotropics are parametrically 
larger than </> 2 . 

The most important 3-point function induced by secondaries on large angular scales is the one 
produced by the correlation between the gravitational lensing effect and the ISW effect. The matter 
distribution along the line of sight both lenses the primary anisotropies and creates additional ones 
through the ISW effect. This results in a non-zero 3-point function first calculated in I12j. Due 
to the enhanced nature of the gravitational lensing effect, this 3-point function is rather large and 
is expected to be detected by the next generation of full sky CMB experiments such as Planck. 

A second source of non-Gaussianity are non-linearities operating at recombination. Not much is 
known about the explicit form of their resulting 3-point function. The naive estimate is that higher 
order corrections are suppressed by powers of the Newtonian potential <j), so that non-Gaussianities 
should be of order 3 • 10~ 5 . It is quite important to refine this estimate. First of all this estimated 
level of non-Gaussianity is roughly of the same order as what the forthcoming Planck mission can 
detect; this implies that to fully exploit Planck data we must be able to control the 3-point function 
generated at recombination. Moreover it is not a priori clear that the naive estimate of a 1CP 5 
non-Gaussianity is a good one; physics at recombination is quite complicated and rich (think about 
the structure of the acoustic peaks in the various spectra) and there could be a sizeable suppression 
or enhancement. 

The calculation of the full 3-point function of the primary anisotropies is very complicated. 
One should take into account all the non-linearities in the evolution of the baryon-photon fluid at 
recombination and non-linearities coming from general relativity. The purpose of this paper is more 

■"■The potential wells are smaller on scales that were inside the horizon during the radiation era as fluctua- 
tions are not able to grow during this period (see for example 
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modest. We are going to derive the result in a particular geometrical limit for the 3-point function. 
In Fourier space we take one of the modes to be of very long wavelength, in such a way that it is 
outside the horizon at recombination. Moreover we take the other two modes to be of much smaller 
wavelength, so that the triangle in Fourier space is very squeezed, with one side much smaller than 
the others. 

There is a simple physical reason why the calculation is much simpler in this limit. At recombi- 
nation the long wavelength mode is outside the horizon so that it cannot change how recombination 
happens; it is completely irrelevant before it reenters in the horizon. After reentering the long 
wavelength mode is physical, so that the 2-point function of the small wavelength modes will be 
different moving along the long mode and this will generate a 3-point function. The background 
mode cannot change the features of the 2-point function which encode physics at recombination (for 
example the ratio between the first and second peek) but it can only change the angle under which 
a given physical scale is observed. Moving along the long wave the 2-point function will be stretched 
or shrunk (generically in an anisotropic way), but will preserve its features. 

In the following section we will study the effect of a long wavelength mode on the 2-point function. 
As the 2-point function is already at first order in the perturbation we have to consider the effect 
of the background wave at first order only. There are different effects but, even though it is simpler 
and more intuitive to study them separately, the separation among them is gauge dependent and 
therefore non-physical. Once we know at leading order in the perturbation the effect of the long- 
wavelength mode, we can calculate the 3-point function for the temperature fluctuations (section |3J). 
The discussion is easily extended to polarization spectra (section EJ) . The significance of our results 
and the comparison with related works are discussed in the conclusions (section EJ) . 



2 The effect of the background wave 

A useful way to write the metric at non-linear order in the presence of scalar perturbations outside 
the horizon is 

ds 2 = a 2 ( V ) \-d V 2 + e^dxidxi] . (1) 

The variable £ is constant outside the horizon |131 ^ because, in the limit in which we can neglect 
gradients, the effect of the exponential is just an irrelevant rescaling of the spatial coordinates. It has 
been shown that the 3-point function of this variable in standard slow-roll inflation is very small, 
being suppressed by a linear combination of the slow-roll parameters: (CCC) ~ C( e ) 7 ?) (CC) 2 - I n the 
following we are going to neglect this original non-linearity as it is much smaller than the effects we 
are going to calculate, which come from the non-linear relation between the temperature fluctuation 
we observe and C- One could prefer to use another non-linear expression for for example the 
variable ( ne w(x) = C( x ) + C( x ) 2 - This new variable would obviously still be conserved outside the 
horizon. But now the 3-point function of Cnew is not suppressed by the slow-roll parameters and it 
would have to be taken into account because of the same order of the effects we are going to calculate. 
Obviously the final physical result, i.e. the 3-point function of the temperature fluctuations, does 
not depend on our choice of variables, but the separation between what is "created" during inflation 
and what derives from the non-linear relation between the metric and what we observe does. 

There is another reason why the exponential form (^Q) is useful. Suppose we want to study, as 
discussed in the Introduction, the effect of a long background mode which is outside the horizon on 
smaller scale fluctuations: ( = £iong + Cshort- Only with the exponential form the background mode 
effect multiplies the short modes and it is thus equivalent to an unobservable coordinate rescaling, 
before the long mode finally reenters in the horizon. This is necessary for our intuitive arguments in 
the Introduction. Taken this useful choice of variables, now we are going to discuss the effect of the 
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background mode. We do it at linear order because it will affect the 2-point function of the short 
modes which is already at first order in the perturbation. 

When the long wavelength mode is still outside the horizon, the metric can be written in terms 
of the Bardeen variable £ as 

ds 2 = a 2 (7]) [-dr] 2 + (1 + 2((x l ))dxidxi] , (2) 

where £ is time independent and slowly varying as a function of x % compared to the Hubble scale. 
This metric is exact up to corrections sub leading in the k/H expansion. As the mode is outside 
the horizon, recombination happens exactly in the same way in different points along the £ wave. 
Moreover, as these coordinates are comoving, recombination happens at the same conformal time 
rj = rjr everywhere 2 . The only effect of the background wave is to rescale the spatial coordinates, so 
that at different points along the £ wave the same physical scale (e.g. the Hubble horizon) will be 
expressed differently in terms of x % : 

a(rj r )(l + C(x l ))Ax = physical scale = const. (3) 

To follow the long wave mode inside the horizon and to compare with the existing literature we 
go to conformal Newtonian gauge when the mode is still outside the horizon. The metric is given 
by 

ds 2 = a 2 (r) [-{I + 2(f){x l ))dT 2 + (I - 2<f){x l ))dx t dx t ] , (4) 

with the well-known relation (ft = — |£ during matter dominance. In the long wavelength limit there 
is no redefinition of the spatial coordinates going from eq. (|2* |) to (| H> . Therefore the x expression of 
a fixed physical distance depends on the position on the background wave as in eq. (j3J). 

In Newtonian gauge there is an additional effect to take into account. Recombination does not 
occur at the same conformal time r everywhere. As in matter dominance a(r) oc r 2 the relation 
between the conformal times in the two gauges is given by 

r=(l-^(x))r ] . (5) 

As recombination happens at a given rj = r/ r everywhere, the expression in Newtonian gauge will be 
given by 

T r (x) = (1 - -<p{x))r] r . (6) 

We are now ready to calculate the effect of the long wavelength mode on the 2-point function. 
Although physics is locally unchanged at recombination by the presence of the long wavelength 
mode, this background wave will change the angular scale at which the 2-point function is observed. 
The final angle of observation is changed by the (ft wave both because the expression in the (r, x) 
coordinates depends on eft and because of the geodesic propagation in the perturbed metric from the 
last scattering surface to us. The separation between these two effects is actually gauge dependent: 
only the total deviation is physical and therefore gauge independent. We will discuss the issue of 
gauge dependence in Appendix A. As we are interested in the angular redefinition at linear order in 
(ft we can consider each effect separately. We start from the (r, x) redefinition. 

Time redefinition. Both time and space redefinition change the angular scale at which a given 
mode is observed. The angular scale is given by 

" AZ (7) 



TO — T r 



2 We are not assuming that recombination is instantaneous. Whatever is the time dependence of the process 
it will be the same everywhere, independently of the position along the long-wave mode. 
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where Ax is a given separation on the last scattering surface and tq — r r is the distance to the last 
scattering surface. The distance to the last scattering surface is dominated by the present conformal 
time To, so that the fluctuations of order (ft in the conformal time of recombination of eq. © change 
the angular scale only by 

66 T r 1 1 

^ (ft^ (ft (8) 

6 T (1 + z r ) 1 ' z 

where z r ~ 1100 is the red-shift at recombination. This effect is therefore much smaller than the 
Ax fluctuations of eq. © which give a relative angular variation of order (ft. Fluctuations in r r can 
therefore be neglected. 

Space redefinition. The long wavelength mode, redefining distances as in eq. gives an 
isotropic (i.e. independent of the orientation with respect to the background mode) redefinition of 
the angle of observation 

T = f*- 

Going to Fourier space, this implies a variation of the Cj's of the form 

8{?C l ) = -\<ft--!^-(?C l ). (10) 
3 a log I 

For a scale invariant spectrum C\ oc l~ 2 there is obviously no effect. 

Let us now consider the effect of the propagation on the perturbed background. There are 
two effects. The first one is the variation of the position of the last scattering surface due to the 
potential along the line of sight (Shapiro time delay), the second is the variation of the angle due to 
the non- homogeneity of the background field (lensing) . 

Shapiro time delay. To study the effect of the Shapiro time delay we take the equation for 
the geodesic in the metric (jlj at zeroth order in the spatial derivatives of (ft, i. e. neglecting lensing 
effects 

r-r r (x) = / (1 -2(ft{x'))dx' . (11) 



This implies that the radial position of the last scattering surface is not only influenced by the 
different time of emission r r (x) of eq. ©, but also by the Newtonian potential <ft(x) along the line 
of sight. This effect has been considered in [T3|. However, the effect is suppressed with respect to 
the fluctuations in Ax discussed above. The reason is that the integral of (ft in the right hand side 
of eq. (|llj) tends to average to zero unless the mode wavevector is perpendicular to the line of sight. 
This implies that this effect will be comparable to eq. © only for the lowest multiples. 

Lensing. A spatially inhomogeneous Newtonian potential (ft will deflect the geodesies and give 
rise to lensing. In the Born approximation the shear matrix $jj can be written as an integral over 
the unperturbed line of sight ^5] 

dSBi n Z^ 3 * D(D lss -D), 



^ = —1 = 2 dD -5^2 '-ViV j( ft(D) , (12) 

OVj Jo ^LSS 

where .Dlss is the comoving distance to the last scattering surface, which we take infinitely thin. 
Although we have a certain cancellation along the line of sight as for the Shapiro time delay, the 
additional gradients acting on the Newtonian potential enhance the lensing effect. In the next 
section we will show that lensing gives a contribution to the 3-point function comparable to the 
space redefinition effect. The part of the shear matrix proportional to the identity matrix gives rise 
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to an isotropic rescaling of the angles, similarly to what happens in eq. ©■ In this case the relative 
variation of the angle from the last scattering surface to us is given by the convergence k 



6e ^lss D(D LSS -D) 

T = K = Jo Dlss ' 

where is the Laplacian restricted to the directions perpendicular to the line of sight. Contribu- 
tions to the shear matrix not proportional to the identity matrix will give an anisotropic redefinition 
of the angle of observation, i.e. dependent on the orientation with respect to the background wave. 

Before going on an important remark. We are going to consider the 3-point function induced by 
lensing, correlating the lensing effect of the background wave with the temperature fluctuations it 
produces close to the last scattering surface. As we mentioned in the introduction, in a Universe with 
a cosmological constant the Newtonian potential is not constant and an additional contribution to 
the temperature fluctuation comes from the Integrated Sachs Wolfe effect. It is well known [ 1 H I12| 
that the correlation of this effect along the line of sight with lensing gives quite a big 3-point function, 
which is sensitive to the value of the cosmological constant. This effect is physically distinct from 
what we are studying as it occurs along the line of sight and it could, at least in principle, be 
separated through an independent reconstruction of the projected potential for example by studying 
weak lensing of galaxies. 



3 The 3-point function 

In this section we explicitly present the 3-point function for the temperature fluctuations coming 
from the effects discussed above, in the flat sky limit. This limit is not a good approximation if the 
background mode has a wavelength comparable to our present horizon, i. e. for its lowest multiples. 
Nevertheless the expressions are much more transparent in this limit so we leave for appendix B the 
full expressions for the spherical geometry. 

We define the spectrum of the 2-point function in the flat sky limit as 

{ iT ^ )ST ^ ) = (^Hk + h)c h . (i4> 

In a similar way we can factor out from the 3-point function the delta of momentum conservation 
defining 

As we discussed in the previous section, the two leading effects of the background wave on the 
2-point function are space redefinition and lensing. The 3-point function from the space redefinition 
effect is obtained correlating the temperature fluctuation given by the long wavelength mode with 
its effect on the 2-point function in eq. ()10j) . We get 

F( r i; r 2; r 3) = <M> £ (-f^eo.)) , m 

where with the prime we indicate the corresponding quantity without the (2ir) 2 5^ of momentum 
conservation. Let us stress again what is the range of validity of this expression. The "background" 
mode l\ must be outside the horizon at recombination and we must be in squeezed, "Maldacena" 
P limit 

h <C 200 , l 2 ^ h > h ■ (17) 
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An expression similar to eq. Q16JI holds for the isotropic lensing contribution 



FffM) = (~fi 2 A)) . (18) 

Note that these two contributions are quite similar and we expect a partial cancellation between the 
two. The reason is that regions close to the last scattering surface with positive convergence n ( 3 ) 
are regions of overdensity, as the Laplacian of <fi is proportional to the mass overdensity. Regions 
of overdensity are regions of negative Newtonian potential 4> and this implies a partial cancellation 
between eq.s (fT^)l and (|T%|) . 

Let us now consider the anisotropic lensing contribution. The background plane wave depends 
only on one direction in the sky, so that in a suitable basis the shear matrix will be of the form 4 

This implies that the anisotropic shear is 

*r = ( S _» ) ■ < 2 °> 

The angle of observation of a given short mode at recombination is redefined in a way that is 
proportional to cos(2i^i2), where (p±2 is the angle between the short mode and long lensing wave. 
This gives us a 3-point function of the form 

F{h;h;h) = K(r i )*) / cos(2^ 12 )i (_^L_(3cy) . (21) 

For any 3-point function we can calculate the signal to noise ratio S/N for an ideal cosmic 
variance limited experiment. It is given in the flat sky approximation by [23] 

For the space redefinition contribution of eq. (fTfj)) the quantity ^ ST ^ ^ 4>{l\)*j can be expressed 

in terms of the power spectrum Pji using the Sachs-Wolfe formula 5T/T = 4>/3, which is valid for 
modes which are outside the horizon at recombination 5 

(M)«5).)'_^ ({ (S)c(Sn'-^^. m 

1/2 

The best constraint on the power spectrum P-jz now comes from the WMAP collaboration P^ ~ 
4.3 x 10 -5 |16j . The logarithmic derivative of the short wavelength spectrum can be calculated 
using publicly available codes like CMBFAST 17 . In figureHwe show the signal to noise ratio as a 

3 In the integral of eq. (jlHJl only the contribution close to the last scattering surface has a substantial 
correlation with the temperature ST(li) of the long wavelength mode. 

4 This is not true for a generic background, like for example a spherical harmonic. 

5 We are neglecting the contribution of the Integrated Sachs- Wolfe effect. This effect is created by the 
variation with time of the potential close to us, at low red-shift. It is therefore correlated to the potential on 
the last scattering surface only for the very low multiples. 
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function of the maximum I used for the short modes in an ideal experiment. We restricted the sum 
to configurations satisfying (|17|) only. 

Let us now concentrate on eq. ()18|) . In the expression for the correlation between the temper- 
ature fluctuation and the convergence of the long wavelength mode we have been sloppy. As the 
convergence is an integral over the line of sight there is no one to one correspondence between a 
mode we see in the sky and the momentum perpendicular to the line of sight. This implies that a 
given mode l\ of 5T/T is not exclusively correlated to a mode — 1\ of the convergence as implicit 
in eq. (|18|) . Nevertheless we have checked that this effect is quite small, except for the lowest mul- 
tiples, as the temperature fluctuation is correlated almost exclusively with the contribution to the 
convergence close to the last scattering surface. In this limit 

/ 6T{ S) ri\*\' 3 f dk z 2tt 2 [ d ™ s dD D D I 



\ T K{1) / 2h Pn .l 2n(k 2 + l 2 /D 2 ss f/ 2 J D LSS D LSS {1 D LSS ] D 2 SS 



;2 



ik z (D-D h33 ) 



(24) 

where we have approximated D = .Dlss is some place. From this we get 

k(T)*)' = ^(~^l 2 )Pn dy —1— J 1 dx x(l - x) cos(y(l - x)) , (25) 

where the integral is easy to evaluate numerically. What is the dependence on I of this expression? 
There is an I 2 enhancement with respect to the correlator IJ23JI coming from the Laplacian in the 
lensing expression. This seems to suggest, neglecting the lowest multiples, that lensing should be 
dominant. Actually the I 2 enhancement is cancelled by two other effects. The first one is the 
cancellation along the line of sight we already discussed for the Shapiro time delay; if the wavevector 
has a component along the line of sight the background wave will lens the 2-point function in opposite 
ways along the path with a resulting partial cancellation. The second effect comes from the window 
function D(Z?lss — D)/D^ss in expression (|T2~)) . This tells us that lensing close to the last scattering 
surface is very ineffective, as it is easy to understand geometrically. The correlation between lensing 
and temperature fluctuation is therefore suppressed. It can be explicitly checked in the integral ()25|) 
that both these effects give an I" 1 suppression, so that the correlator (|25fl finally behaves like the 
one in eq. (|23j) . it goes as l~ 2 . Space redefinition and lensing give comparable contributions. 

Expression (|25j) allows us to estimate S/N for the 3-point function of eq. (|18[) . In figure Q we 
show the result as a function of the maximum I for an ideal experiment. We also show S/N for the 
sum of the two effects discussed so far to show that there is a partial cancellation between the two. 

For the anisotropic lensing, it is easy to relate its signal to noise ratio to the one given by the 
isotropic lensing effect. Averaging expression (|21[) over the orientation between the background mode 
and the 2-point function we get a signal to noise ratio smaller by a factor of \/2 with respect to the 
isotropic lensing effect. We plot the result in fig. H As a consequence of the partial cancellation 
between the isotropic effects, the anisotropic lensing gives the biggest signal to noise. 



4 Polarization 3-point functions. 

We can easily extend the results of the previous sections to the 3-point functions involving polariza- 
tion modes. Note first of all that the E polarization decays quite fast for modes outside the horizon 
at recombination, while 5-modes are zero in absence of a tensor contribution. This implies that, in 
the limit we are considering {i.e. one long- wavelength mode outside the horizon at recombination 
and two much shorter modes), the long mode must be a temperature mode. Therefore we have only 
to consider four 3-point functions involving polarization modes: TTE, TEE, TTB, TEB, where 
the first mode is the long-wavelength mode. 
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Figure 1: Signal to noise as a function of the maximum I for an ideal, cosmic variance limited 
experiment. Dashed line: space redefinition effect. Continuous thin line: isotropic lensing effect. 
Thick line: total effect. We see that the total signal to noise is smaller than the individual contri- 
butions as the two effects partially cancel. The parameters of the cosmological model we used were: 
n b = 0.045, n c = 0.255, ft A = 0.7, h = 0.7 and a 8 = 0.85. 
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Figure 2: Signal to noise for the anisotropic lensing contribution (thick line) compared to the sum 
of the isotropic contributions (thin line) as a function of the maximum / for an ideal, cosmic variance 
limited experiment. The anisotropic effect has a bigger signal to noise. Cosmological parameters 
are the same as in the previous figures. 



The discussion for the TTE and TEE modes follows closely what we said about the 3-point 
function of the temperature map. The effect of the background mode is always to redefine the angle 
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at which a certain physical scale is observed. Given the partial cancellation discussed above of the 
two isotropic effects (space redefinition and isotropic lensing) the leading effect is the anisotropic 
lensing. The explicit form for the TTE and TEE correlators is the same as in eq. ([21 j). with the 
only difference that the small scale spectrum Ci 2 is respectively the TE and the EE spectrum: 



where X stands for either T or E. 

It is straightforward to calculate the signal to noise of these correlators from eq. 



(26) 



(27) 



(q/atV 1 f d2l id 2 l 2 F TTE (h;l 2 ;l 3 ) 2 2 _ I f d 2 hd% F TEE (h;l 2 ;h) 2 

WVTTE-~J (2yr)2 2C TT C TT C EE W»)TEE-~J ^2 ^qTTqEEqEE 

In calculating the variances at denominator in the first expression we neglect contributions involving 
(Cf E ) 2 because it is much smaller than Cf T ■ C, , i.e. the correlation between E and T is not 
perfect and it becomes quite small at large I. 

The results are given in figure Has a function of the maximum I for an ideal cosmic variance 
limited experiment. 




1000 1500 2000 2500 3000 

Figure 3: Signal to noise for the TTE (continuous line) and TEE (dashed line) correlators as a 
function of the maximum I for an ideal, cosmic variance limited experiment. Cosmological parameters 
are the same as in the previous figures. 



We can finally discuss correlators involving B modes. Even if we do not have any original B 
mode from a tensor component, the anisotropic lensing generates a B component from a pure E 
polarization. Lensing does not change the state of polarization of a given photon, but changes the 
angle at which we observe a given physical scale. The magnitude and orientation of the wavevector of 
a given mode is therefore changed. A pure E mode is converted into a mixing of the two polarization 
states. Taking an E mode in the flat sky limit it is straightforward to see that the amount of B 
polarization generated by anisotropic lensing is given by 

B = -EKsm(2cp 12 ) , (28) 



10 



where k is the convergence and <pi2 is the angle between the lensing long-wavelength mode and the 
short E mode. 

Prom this we get the expression for the TTB and TEB correlators 



Fttb{Ii-M-M) = ~( 6 -^ K(hy)' sm(2 Vl2 )C? 2 E , 
Fteb{Ii-M-M) = -2(^^ K (r i )*) / sin(2 9 9 12 )C^. 



(29) 



(30) 

The factor of 2 difference between the two expressions comes from the two possible E to B conversions 
of the EE correlator. Note that here we do not have any derivative of the short wavelength spectrum: 
the effect of the long mode is to create a B polarization and not to just stretch or shrink a preexisting 
2-point function. 

The signal to noise ratios are 



{S/Nf TTB 



1 



7T 



d 2 hd 2 l 2 F TTB (h;r 2 ; 



2Cl T Cl T C* B 



(S/N) 



TEB 



1 



7T 



<Ph<PhF TEB (h;h;h? 
(2*)' CjTC?*C?s 



(31) 



For the variance Cf B we take the one created by the lensing itself. The smallness of the BB 2-point 
function gives a bigger signal to noise to these correlators involving B modes with respect to the 
others, at least in the case of an ideal experiment, limited only by cosmic variance. Moreover as the 
BB 2-point function coming from lensing peaks at I ~ 1000, the S/N rises fast and then saturates 
around I ~ 1500. 
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Figure 4: Signal to noise for the TTB (continuous line) and TEB (dashed line) correlators as a 
function of the maximum I for an ideal, cosmic variance limited experiment. Cosmological parameters 
are the same as in the previous figures. 



5 Conclusions 

In this paper we calculated the 3-point function for the CMB temperature and polarization maps 
generated by non-linearities at recombination, in the limit in which one of the modes is outside the 
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horizon at recombination and the other two are of much smaller length. These generated 3-point 
functions, besides their intrinsic interest, represent some sort of "background" if we are interested 
in studying the non-Gaussianity of the primordial spectrum of perturbations. Even though dif- 
ferent models predict different shape dependences for the primordial 3-point function ^H], usually 
experimental limits are set on the scalar variable /nl defined through the relation 

C(x) = Q(x) - lfNL(C g (x) 2 - (Q) , (32) 

where £ is the observed perturbation and Q g is gaussian ( 6 ). The present limits on /nl are given by 
the WMAP experiment —58 < /nl < 134 at 95% of CL |19j . The Planck experiment will improve 
this limit down to /nl ~ 5. 

A useful way to express our result on the signal to noise of the TTT 3-point function is to 
calculate the equivalent /nl which gives the same signal to noise. Summing only over the triangles 
for which our approximation applies up to l max = 3000 we get 

/Nf" ~ 0.7 . (33) 

This means that this non-Gaussianity should not be detected by Planck. Any signal of 3-point func- 
tion (besides the well known ISW/lensing correlation) should be considered of primordial origin, 
at least in the geometrical limit in which our calculation applies. It is possible that future experi- 
ments will detect the 3-point functions generated by non-linearities at recombination especially using 
polarization modes, for which we obtained a bigger signal to noise ratio. 

We stress that the configuration dependence of our result is completely different from that 
implied by eq. (|32j) which, in the geometrical limit we are considering, is proportional to the 2-point 
spectrum of the short modes. Our 3-point function, on the other hand, oscillates around zero as 
a function of the length of the short modes, being essentially proportional to the derivative of the 
2-point spectrum. That is to say the non-Gaussian signature produced by the /nl model in the 
presence of a long wavelength mode is a change in the amplitude of the 2-point function on small 
scales, a shift up and down of the Ci curve. The effect we studied here changes the angular scales of 
the features in the 2-point function, a sideways shift in the C\ curve. Moreover, as the anisotropic 
effect is dominant, the 3-point function depends on the relative orientation between the long and 
short modes. These differences should further reduce the contamination if we are looking for a 
primordial /nl- 

Even though our calculation was done in a particular geometrical limit for the modes entering in 
the 3-point function, the signal to noise we get summing over these particular squeezed configurations 
is not very small with respect to what is expected for the full 3-point function. The reason is that, 
in the limit of a scale invariant 2 and 3-point function, the signal to noise of eq. (|22|) diverges by 
dimensional analysis as Z^ ax (in eq. (|22j) F oc and C oc l~ 2 ), with a possible additional log 
divergence which depends on the particular form of the 3-point function. Even restricting to the 
limit in which our results apply S/N goes as Z^ ax so that we expect to lose only a logarithmic factor 
by our inability to do the full calculation. This estimate does not consider factors of order one given 
by the physics at recombination which could enhance or suppress substantially S/N in the regions 
in which our approximation is not valid. 

References |201 calculate the TTT 3-point function in the limit in which all modes are outside 
the horizon at recombination. The result should be compatible with ours in the limit of squeezed 
triangles, with one side much bigger that the others, with all the modes outside the horizon. The 

6 The non-linearity parameter is defined for the Newtonian potential in matter dominance; this explains 
the factor of 3/5 in the previous expression. 
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3-point function should go to zero because the spectrum is flat at large scales so that its derivative 
entering in eqs (|18|) and (|21l) vanishes. The result of ref.s [2UI 121] does not show this behaviour. 
We were unable to trace back the origin of the disagreement between the results as the two approaches 
to the problem are completely different. 
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The 3-point function is obtained in the main text starting from the variation of the 2-point function 
moving along the background mode. The reader could be worried about this procedure because it 
assumes that we can define the 2-point function before the long wavelength mode reenters in the 
horizon. The usual derivation of the Sachs- Wolfe effect involves an integral over the photon trajectory 
from the last scattering surface to us. This seems to indicate that the temperature fluctuation cannot 
be defined until the present time. 

This is not quite true. Imagine that we have friends in different directions who look outwards 
along our line of sight. Imagine that they are sufficiently close to our last scattering surface to 
observe the 2-point function of the short modes on a small patch of it, just before the background 
mode reenters. We should specify which is the state of motion of our friends or, equivalently, the 
coordinate system with respect to which they are stationary. In the unperturbed Universe we could 
take them to be stationary with respect to the cosmic fluid, but including perturbations, different 
gauge choices will give a different relative velocity of order <p to our set of friends. This will only 
change at leading order the dipole they observe, but if we are interested in higher multiples the 
result is unique and gauge independent. All this set of friends will observe, in their different patch in 
the sky, a 2-point function with the same statistical properties because the background mode is still 
outside the horizon and therefore irrelevant: it just induces an unobservable rescaling of the spatial 
coordinates between one friend and another. From this point of view the effect of the background 
mode is to induce a deformation between the map one of our friends observed and the map of the 
same patch on the last scattering surface we observe now. The effect is physical and gauge invariant. 

Thus the full effect is gauge invariant, but as we stressed in the main text, the separation of 
the various effects we discussed is not physical and depends on the gauge we choose. To show this 
in a concrete example let us start from the conformal Newtonian gauge (@J and redefine the time 
slicing to make it comoving, without touching the space coordinates (total matter gauge) . In matter 
dominance the time shift required is given by 



The transformed metric is given as a function of the new time variable f by 



Appendix A. Gauge invariance 




(35) 



Outside the horizon, neglecting the gradient of <f> and using the relation 4> = — 1£, we get the metric 
© with f = r]. 
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As the gauge we are considering is comoving, there is no time redefinition effect in this metric; 
recombination happens at the same time f r everywhere. As there is no redefinition of spatial coor- 
dinates, the space redefinition and lensing effects are obviously left untouched. Therefore the time 
redefinition effect must be incorporated in the Shapiro effect in this gauge. It is easy to verify this. 
At leading order in <j) and neglecting lensing the null geodesies in the metric (|35j) follow the equation 

df = (l --0 + -TV(j))dx . (36) 
3 3 

This gives integrating by parts 

f - f r = J (1 - 2(f>(x , ))dx' - f r ^{x) , (37) 

which coincides with equation The fluctuation along the background wave of the time of 

recombination in a generic gauge becomes a contribution to the distance to the last scattering 
surface in a comoving gauge. 

The lensing effect is, on the other hand, independent of the parameterization of the perturbations. 
Propagating backwards from us the null geodesies, we can covariantly define the shear tensor on a 
surface perpendicular to the line of sight as the pull back on the surface of the covariant derivative 
[22] 

Vpk„ , (38) 

where is the null vector field tangent to the geodesies. As we are interested in the result at 
first order in the perturbation, both the surface perpendicular to the line of sight and the affine 
parameterization of the geodesies can be made in the unperturbed Universe and it is therefore 
common to all gauges. The integral of the shear tensor along the line of sight is therefore gauge 
independent. 



Appendix B. Expressions on the full sky 

In the main text we discussed our results in the flat sky limit, which is much more transparent. In 
this appendix we generalize some of the expressions to full sky. The correspondence between flat 
sky and full sky for the 3-point functions is thoroughly studied in |23j . 

In a full sky analysis definition Q15JI . which in the flat sky limit enforces the translational invari- 
ance of the 3-point function, is replaced by 

((6T/T) hmi (ST/T) l2m2 (5T/T) hm3 ) = ( ^ ^ ^ \ B(l h l 2 ,l 3 ) , (39) 

involving the Wigner 3j symbols, which is manifestly rotationally invariant. 

It is straightforward to extend the flat sky 3-point functions to the full sky. In the case of 
isotropic contributions, in which there in no dependence on the relative angle between the long and 
the short modes 

*.«.,«.,« - ( To o ) ^ + 1)( ^ 1)( * + V .;S,S) , (40) 

where F in eq.s (|16|) and (|18j) must be written as product of the two full sky spectra. The previous 
expression is different from zero only for l\ + 1% + Z3 even, as implied by parity invariance. 
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For the anisotropic contribution eq. (|21|) the cosine of the angle between the long and short 
modes goes into the Wigner 3j symbols giving 



o n j n _ h h h \ (2/ 1 + l)(2/ 2 + l)(2/3 + l) /JT \ 1 / 3 2 A 

B^ {hM,h) -[ 2 _ 2 J V ^ \ T Ail TaloiT^V • 

(41) 

This expression holds only for l\ + I2 + h even, while it is taken to be zero for parity violating 
correlations. The same correspondence holds for the polarization spectra TTE and TEE eq. (|26|). 
For the 3-point functions involving B modes we have 

Brr^MM) = i ( \ \ % ) y (2 '' + 1)( \ +1)(2 ' 3 + 1) (f C,r (42) 

and similarly for the TEB correlator. In this case the previous expression holds for l\ + I2 + I3 odd, 
while it is taken to be zero in the other cases. 

In the full sky the integral over the line of sight which gives the temperature-convergence corre- 
lation (see eq. (|24|) for the flat sky expression) takes the form 

((5T/T) hmi n\ m \ = 8 mim2 5 hh ■ ^(-4vr)/(/ + l)P n H -jlM [ — (1 " x)j h {xy) , (43) 

N ' ^0 jo y Jo x 

where j are spherical Bessel functions. In the spherical geometry we do not have the problem we 
discussed in the main text about the correspondence between the transverse momentum and the 
mode we observe in the sky. We can check that the flat sky approximation used in the main text is 
quite good, except for the lowest multiples. 
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